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The FAA has worked with Transport Canada and others to develop allowance times for 
aircraft operating in ice-pellet precipitation. Wind-tunnel testing has been carried out to 
better understand the flowoff characteristics and resulting aerodynamic effects of anti-icing 
fluids contaminated with ice pellets using a thin, high-performance wing section at the 
National Research Council of Canada Propulsion and Icing Wind Tunnel. The objective of 
this paper is to characterize the aerodynamic behavior of this wing section in order to better 
understand the adverse aerodynamic effects of anti-icing fluids and ice-pellet contamination. 
Aerodynamic performance data, boundary-layer surveys and flow visualization were 
conducted at a Reynolds number of approximately 6.0xl0 6 and a Mach number of 0.12. The 
clean, baseline model exhibited leading-edge stall characteristics including a leading-edge 
laminar separation bubble and minimal or no separation on the trailing edge of the main 
element or flap. These results were consistent with expected 2-D aerodynamics and showed 
no anomalies that could adversely affect the evaluation of anti-icing fluids and ice-pellet 
contamination on the wing. Tests conducted with roughness and leading-edge flow 
disturbances helped to explain the aerodynamic impact of the anti-icing fluids and 
contamination. The stalling characteristics of the wing section with fluid and contamination 
appear to be driven at least partially by the effects of a secondary wave of fluid that forms 
near the leading edge as the wing is rotated in the simulated takeoff profile. These results 
have provided a much more complete understanding of the adverse aerodynamic effects of 
anti-icing fluids and ice-pellet contamination on this wing section. This is important since 
these results are used, in part, to develop the ice-pellet allowance times that are applicable to 
many different airplanes. 


Nomenclature 


b 

= model span 

c 

= chord length 

Cd 

= wing drag coefficient 

C l 

= wing lift coefficient 

Cm 

= wing quarter-chord pitching moment 

Cl, max 

= maximum lift coefficient 

H 

= boundary-layer shape factor 

k 

= roughness or disturbance height 
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I. Introduction 

The proper use of deicing and anti-icing fluids is required to maintain adequate safety margins for takeoffs in 
cold-weather-precipitation conditions. The FAA and other civil airworthiness authorities publish holdover time 
tables (HOTs) that determine, in part, the length of time airplanes may remain on the ground after treatment with 
deicing and anti-icing fluids depending upon the weather conditions. HOTs are determined based on the results of 
endurance time testing of ground fluids in natural or simulated icing conditions. The endurance times are 
determined by applying failure criteria using visual cues, which vary depending on the weather conditions (snow, 
freezing drizzle, freezing rain). While HOTs exist for most types of weather conditions, ice-pellet precipitation and 
ice-pellets mixed with other types of precipitation 1 have demonstrated different behaviors. The difficulty in 
establishing visual cues for detecting fluid failure in ice-pellet conditions has thus far precluded the establishment of 
practical failure criteria for endurance time testing. The FAA has never provided HOTs for operations in ice-pellet 
conditions since HOTs cannot be established without endurance time testing. 

In 2005, the FAA became concerned about operations in ice-pellet conditions after learning of questions raised 
relating to regional carriers operating in ice pellets with non-slatted aircraft and with less experienced flight crews 
than the major carriers. As a result of these concerns, the FAA issued three notices effectively disallowing 
operations in ice pellets in October 2005. 2-4 In the winter of 2005-2006, this action had a major impact at several 
airports during ice-pellet events. Due in part to concerns raised by these events, the FAA undertook a thorough 
reassessment of its policy with respect to operation in ice-pellet conditions. Environmental chamber testing also 
raised questions regarding the ability of the anti-icing fluid to flow off aerodynamic surfaces during a takeoff run. 
Adhering contamination or a significant amount of fluid remaining on the wing at takeoff would both present a 
safety hazard. 

The significant differences in ice-pellet precipitation behavior relative to other forms of precipitation combined 
with the potential flight safety hazards led the FAA Ground Icing Program and Transport Canada to sponsor aircraft- 
based and wind-tunnel based research to better determine allowance times associated with operations in ice-pellet 
conditions. An allowance time differs from a HOT in that there is no provision for a pre-takeoff contamination 
check if the time is exceeded, meaning that the allowance time cannot be extended by visual or external tactile 
checks of the aircraft critical surfaces within five minutes of takeoff. Ordinarily pre-takeoff contamination checks 
are conducted from the cabin, and the research showed that such a check would be of little use in the case of ice 
pellets. 5 Broeren and Riley 6,7 provide a more detailed summary of this flight-based and wind-tunnel based research 
along with a description of the resulting ice-pellet allowance time tables. 8 " 10 Of interest here is the evolution of the 
wind-tunnel experiments and resulting safety implications. 

A series of test campaigns were conducted beginning in January 2007 through February 2012 at the Propulsion 
and Icing Wind Tunnel (PIWT) owned and operated by the National Research Council (NRC) of Canada. 11,12 Ice 
pellets were simulated using crushed ice from blenders which was run through sieves yielding a size distribution 
similar to natural ice pellets. An airplane takeoff acceleration profile and rotation were simulated in the wind-tunnel 
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testing. Initially, allowance times were determined using the visual observations of the researchers and existing 
single-element airfoil models were used for testing. Beginning in January-February 2010, tests were continued in 
the NRC PIWT using a thin, high-performance wing section rigged with a single-element slotted flap fixed at 20 
deg. The thin, high-performance airfoil was regarded as typical of wings on regional jet transport airplanes. 
Aerodynamic measurements were again made, and effects on lift coefficient were significantly larger than in past 
years. Along with visual evaluation by the researchers, a criterion for percentage loss in lift at 8 deg. angle of attack 
was explicitly used in developing allowance times. Based on the results, the FAA published an updated ice-pellet 
allowance time table for winter 2010-201 1 and 201 1-20 12. 13 

The introduction of an aerodynamic-performance based criterion and subsequent changes to the ice-pellet 
allowance time tables led to concerns regarding the fidelity and applicability of the measured lift degradation in the 
ice-pellet contamination fluid testing at the NRC PIWT. Given the potential operational safety implications 
associated with the measured performance degradations, it is even more important that these concerns be addressed 
appropriately. The essence of these concerns is understanding how the performance degradations measured on the 
two-dimensional model in the NRC PIWT testing relate to an actual three-dimensional, full-scale airplane 
configuration. Broeren and Riley 6,7 developed a scaling method based upon the Aerodynamic Acceptance Test 
(AAT) that scaled the percent lift loss on the thin-high performance wing to the percent loss in maximum lift on a 
full-scale Boeing 737-200ADV airplane. This scaling relationship allowed the aerodynamic performance penalty to 
be used as one of several criteria in determining allowance times for ice-pellet precipitation conditions. 

Broeren and Riley 6,7 also addressed a number of concerns regarding the aerodynamic testing conducted at PIWT. 
These issues were related to the design and conduct of the experiments relative to standard aerodynamic test 
methodology. While this review showed that the tests were conducted in accordance with commonly accepted 
practices, no baseline aerodynamic characterization of the thin, high-performance wing section had been conducted. 
Due to the unique design of the model and lack of any comparable data, such a characterization was recommended 
in order to better understand the airfoil aerodynamic behavior. Research of this type had been carried out by Hill 
and Zierten et al. 14-17 in their development of the AAT that is used to qualify deicing and anti-icing fluids. 18 The 
authors used both 2-D and 3-D models and conducted extensive flow visualization studies on the clean, dry models 
as well as aerodynamic performance studies with roughness and simulated frost applied to the models. Data from 
the latter experiments were used to help explain the adverse performance effects observed during the fluid tests and 
relate these adverse effects to those observed on larger-scale models. A similar approach was followed in the 
present research using the thin, high-performance wing. The present research was necessary because the design of 
the thin, high-performance wing is significantly different from airfoils and wings used in past research. 
Furthermore, previous work was conducted before the advent of Type IV anti-icing fluids used in the present 
research along with the additional effects of ice-pellet contamination. These factors required additional research 
into the aerodynamic effects of the fluids and contamination. 

The objective of the present research was to characterize the aerodynamic behavior of the thin, high-performance 
wing used for aerodynamic testing of ground anti-icing fluids. Aerodynamic performance measurements, flow 
visualization and boundary-layer surveys were conducted on the clean, dry wing. Tests were performed with 
leading-edge roughness (simulating ice roughness) and roughness applied to the entire upper surface (simulating 
frost). Additional tests were conducted to simulate the specific effects of the transitory nature of the anti-icing 
fluids. This characterization is important to understand the adverse aerodynamic effects of anti-icing fluids and ice- 
pellet contamination on this model since the results are used, in part, to develop ice-pellet allowance times that are 
applicable to many different airplanes. 


II. Experimental Methods 

All of the experiments were carried out using the thin, high-performance wing section in the NRC PIWT facility 
located on the Montreal Road campus in Ottawa. The PIWT is an open-return wind tunnel that draws air from the 
outdoors and is therefore naturally cooled. The fan is located in the wind-tunnel inlet and was driven by a gas 
turbine engine for the test campaign. The test section measures 10 ft (3 m) across and 20 ft (6 m) in height. An 
insert was used that reduced the height to 16.4 ft (5 m) resulting in an increase in maximum test-section speed. The 
wing was mounted horizontally in the test section as shown in Fig. 1. The model had a span of 7.9 ft (2.4 m) and 
was isolated from the wind-tunnel walls with endplates. There were large gaps between the endplates and the test- 
section walls (~ 1 ft (0.3 m) on each side). The size of the endplates was minimized to facilitate visual access to the 
wing surface which was very important for the experiments with anti-icing fluids and contamination. An additional 
motivation of this work was to determine any adverse effect of these gaps coupled with the relatively small size of 
the endplates. Figure 2 shows the 6 ft (1.8 m) chord, thin, high-performance wing section with a single-slotted flap 
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deflected at 20 deg. The model had a fixed (or “hard”) leading edge and the airfoil design was regarded as 
representative of wings on regional jet transport airplanes. 

The model was supported by two-component external force balances at each end, located outside of the test- 
section walls. Each balance rotated with the wing section, and consisted of two load cells in the normal direction 
and one load cell in the axial direction. A drive motor and gearbox connected to the left-side force balance provided 
rotation for setting the wing geometric pitch angle. Lift, drag and pitching moment data were measured from the 
balance during angle of attack sweeps. These data were corrected for wind-tunnel wall effects using the 2-D airfoil 
methods described in Rae and Pope. 19 Performance data were acquired over a small range of Reynolds and Mach 
number depending upon the outside air temperature and the speed. These conditions are summarized in Table 1. 
This experimental set-up is very similar to that used to evaluate the performance of fluids in their uncontaminated 
state and with ice-pellet contamination and has been described elsewhere. 20,21 



Figure 1. Photograph of thin, high-performance wing section installed in the NRC PIWT facility and end- 
view drawing of endplate configuration. 


0.10 i 

0.05 ~ 



Figure 2. Thin, high-performance airfoil section used in the PIWT testing, reference chord = 6.0 ft (1.82 m). 


Table 1. Reynolds and Mach Number Conditions 


Speed, knots 

Temperature, °C 

Reynolds Number 

Mach Number 

80 

0 

5.6xl0 6 

0.12 

80 

-20 

6.5xl0 6 

0.13 

100 

0 

7.0xl0 6 

0.16 

100 

-20 

8.1xl0 6 

0.16 
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A detailed uncertainty analysis was performed for the aerodynamic data for 20:1 odds. The combined non- 
linearity and repeatability for each load cell in the balance is ±0.11% of the rated output. This results in an 
uncertainty of ±4.47 lbf (±19.9 N) in the total normal force and ±1.25 lbf (±5.6 N) in the total axial force. The 
uncertainty associated with the pitching moment about the balance center is ±6.13 lbf- ft (±8.31 Nm). The dynamic 
pressure in the test section was calculated based on a pre-defmed calibration constant and the measured dynamic 
pressure in the wind-tunnel contraction. The contraction dynamic pressure has a combined precision and accuracy of 
±4.38><10" 5 psi (±0.302 Pa) that corresponds to an uncertainty in the test section dynamic pressure of ±0.0419% at 80 
knots and ±0.0283% at 100 knots. The combined uncertainties of the balance loads and the test-section dynamic 
pressure were used to calculate the instrumentation error associated with each of the aerodynamic coefficients as a 
function of the test-section velocity and the wing geometric pitch angle. The uncertainty associated with the 
geometric pitch angle is based on the accuracy of the inclinometer used to set the wing pitch angle and the backlash 
in the gear box system that drives the wing rotation. In addition to instrumentation error, the repeatability of each 
measurement point was considered in the evaluation of the overall data uncertainty. For the clean-wing cases, the 
standard deviation of the aerodynamic coefficients at each angle of attack was calculated from numerous repeat runs 
conducted throughout the test campaign. The instrumentation errors and the standard deviations of the clean wing 
measurements were combined to provide the total uncertainties associated with the lift, drag and pitching moment 
coefficients. The average of the total uncertainties for -2° < a < 21° are shown in Table 2. All of these uncertainties 
were considered acceptable given the objectives of this research. 


Table 2. Summary of Aerodynamic Performance Uncertainty 


Quantity 

Total Uncertainty 
at 80 knots 

Total Uncertainty 
at 100 knots 

Cl 

±0.012 

±0.011 

C d 

±0.004 

±0.004 

Cm 

±0.002 

±0.002 


Aerodynamic performance data from the force balance were supported by flow visualization and boundary-layer 
measurements. Flow visualization was conducted using thread tufts and surface oil in separate experiments. Tufts 
were applied across the span of the model in several chordwise rows on the main element and flap. The tufts 
provided some information about the unsteady characteristics of the flow as the airfoil model was pitched through 
stall. The surface oil flow visualization provided a higher-resolution, but time-averaged, view of the boundary-layer 
state on the model upper surface. Standard grade mineral oil was dyed with black carbon powder and applied to the 
model with a paint roller. The bare aluminum surface of the model was covered with a smooth, self-adhesive, 
plastic film to provide a uniform background. After the oil was applied, the angle of attack and tunnel speed were 
set to the desired values. Steady state was achieved after a few minutes of run time and was confirmed by visual 
observation through the test-section windows. Photographs were taken at the end of each run. 

The boundary-layer velocity profiles were measured using a rake with 21 total pressure tubes as shown in Fig. 3. 
This rake also contained a single tube used to measure the local static pressure that was used to compute the mean 
streamwise velocity from the total pressures. The rake was placed at three spanwise locations near the trailing edge 
of the main element. These locations were near the midspan station and approximately 2 ft. (0.6 m) to the right and 
left of the midspan station (y/b = ±0.25). Boundary-layer profiles were also measured near the trailing edge of the 
flap at the midspan station. The mean streamwise velocity profiles were subsequently integrated to obtain the 
boundary-layer displacement thickness, momentum thickness and shape factor using the following equations for 
incompressible flow: 





[Eq. 1] 



u 

U edge 



-*-) 

U edge/ 


dz' 


[Eq. 2] 


H = j [Eq. 3] 

The boundary-layer pressures were measured using a ZOC 33 module that had a 1-psid (6,895 Pa) range and an 
accuracy of ±0.001 psi (±6.9 Pa). The resulting uncertainties in typical boundary-layer velocities ranged from 
±1.21% for those closest to the wing surface to ±0.32% for those at the boundary-layer edge — assuming a 
freestream speed of 100 knots. 
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Figure 3. Pressure rake used to measure boundary-layer velocity profiles shown at midspan station near 
trailing edge of main element. 


A number of experiments were conducted with simulated leading-edge ice roughness and upper-surface frost. 
These sources of environmental contamination were simulated with sandpaper roughness of various grit size as 
listed in Table 3. Three roughness sizes representing initial ice accretion were applied to the leading edge from x/c = 
0.08 on the lower surface to x/c = 0.05 on the upper surface. The k/c = 0.00010 roughness representing frost was 
applied over the entire upper surface starting at x/c = 0.08 on the lower surface. Following the method in Runyan et 
al., 17 the k/c = 0.00010 roughness was systematically removed from the leading edge toward the trailing edge. This 
provided an approximate simulation of the aerodynamic effect of the anti-icing fluid as it flowed off of the airfoil 
model. Other tests were performed with roughness applied to the flap only. The effect of leading-edge disturbances 
was also investigated using layers of aluminum speed tape applied across the span of the model near the leading 
edge. Table 4 summarizes the disturbance heights that were investigated. 


Table 3. Sandpaper Roughness Sizes 


Grit Size 

Roughness Height, k , in. (mm) 

k/c 

40 

0.0166 (0.42) 

0.00023 

80 

0.0072 (0.18) 

0.00010 

150 

0.0036 (0.09) 

0.00005 

Table 4. Leading-Edge Disturbance Sizes 

Tape 

Tape 


Height, k (in) 

Height, k (mm) 

k/c 

0.020 

0.51 

0.00028 

0.030 

0.76 

0.00042 

0.040 

1.02 

0.00056 

0.050 

1.27 

0.00069 


III. Results and Discussion 

A. Clean Model Aerodynamics 

The aerodynamic characteristics of the clean wing model were measured to confirm the airfoil stall type, show 
that no flow anomalies existed and document the state of the boundary layer near the airfoil trailing edge. For 
typical anti-icing fluids testing, procedures were established, and implemented in this research, to simulate an 
airplane takeoff profile. The angle of attack was set to -2 deg. before accelerating the test section airflow from 20 
knots (10.3 m/s) to 100 knots (51.4 m/s) in less than 31 seconds. In a typical profile, the model was rotated from a = 
-2 deg. to 8 deg. at a rotation rate of 2.7 deg./sec. after the air speed reached 100 knots. The model was held at 8 
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deg. for 1.5 seconds to gather more data at this angle of attack, then the angle was reduced. In other runs, the angle 
of attack was increased through stall along with other variations in rotation speed and time. This procedure required 
conditional (or phase) averaging of the data during post-processing to develop lift, drag and pitching moment 
information versus angle of attack. These results were compared to aerodynamic performance data that were 
acquired at fixed angle of attack in a pitch-pause mode with one-degree increments in angle of attack. As shown in 
Fig. 4, there was some difference in the maximum lift characteristics depending upon the angle of attack sweep rate. 
The actual value of C L>max varied only from 2.20 to 2.24 while a sta u varied from 19.5 to 20.5 deg. These differences 
may be explained in part by unsteady aerodynamic effects. The conditional averaging procedure required for the 
acquisition of dynamic data can also be affected by the magnitude of the unsteadiness versus the number of samples 
available to be averaged. Therefore some variation in the dynamic data is expected. Repeat runs in pitch-pause 
mode yielded a higher degree of repeatability in the performance data, particularly at stall. Therefore, all of the data 
acquired for these aerodynamic characterization experiments were carried out using the pitch-pause mode. 


2.4 
2.2 
2.0 
1.8 

C L 

1.6 

1.4 
1.2 
1.0 
0.8 

0.6 ■ i ■ i ■ i ■ i ■ i ■ i ■ i ■ i ■ i ■ i ■ i ■ i ■ i ■ kj.\j • i • i • i • i ■ i • i • i • i ■ i ■ i ■ i • i ■ i ■ i 

-4 -2 0 2 4 6 8 10 12 14 16 18 20 22 24 -4 -2 0 2 4 6 8 10 12 14 16 18 20 22 24 

a (deg) a (deg) 



Figure 4. Effect of angle of attack sweep rate on airfoil aerodynamic performance. 


One of the characteristics of the thin, high-performance wing section exhibited in Fig. 4 is the sharp decrease in 
lift coefficient at stall. This behavior is typical of leading-edge stall which has been defined as abrupt flow 
separation near the leading edge without subsequent reattachment. 22 The “abrupt” separation usually results from a 
small leading-edge laminar separation bubble which “bursts” at stall and causes the sharp decrease in lift. 22 Flow 
visualization was conducted to confirm the leading-edge stall type of the airfoil and to determine if there were any 
adverse effects due to the unique model installation with the relatively small endplates. Surface oil flow 
visualization indicated that there was a small, leading-edge laminar separation bubble located on the upper surface. 
This situation is depicted in Fig. 5 with the model at 8 deg. angle of attack. Figure 5 also shows a phenomenological 
sketch of the flow that was adapted from Tanner 23 but has been documented in numerous other references (e.g., 
Ref. 24-27). The separated flow region inside the bubble is clearly visible in the photograph as the thick black line 
of accumulated surface oil. The boundary-layer reattachment region was located immediately downstream at 
approximately x/c = 0.05. Increasing the angle of attack to 18 deg., which is about two degrees prior to stall, had the 
effect of moving the laminar separation bubble upstream as depicted in Fig. 6. Although very difficult to judge from 
this photograph the reattachment region is estimated to be at approximately x/c = 0.03. An additional observation 
from the flow visualization images in Fig. 5 and 6 is that the separation bubble had a very clear two-dimensional 
structure across the span of the model which indicates that two-dimensional flow was maintained on the leading- 
edge area approaching aerodynamic stall. 
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x/c = 0.05 

- Oil accumulation region 


Flow 



Figure 5. Photograph of surface-oil flow visualization (top) at a = 8 deg. and phenomenological sketch of 
laminar separation bubble flowfield, adapted from Tanner. 23 
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Verifying the extent of two-dimensional flow approaching stall was an important part of this work because the 
influence of the small endplates and large gaps was unknown. Initially, it was expected that the combination of the 
small endplate size and relatively large gaps between the test-section walls and the model endplates would lead to 
significant three-dimensional effects near stall. Additional flow visualization images were acquired for the other 
portions of the wing model at a = 18 deg. which was determined to be sufficiently close to the stalling angle to 
characterize the surface flow. Shown in Fig. 7 are surface-oil flow visualization images near the trailing edge of the 
wing main element extending to the left and right-side endplates. The images show that the boundary-layer 
remained attached and that the surface flow was highly two-dimensional, even in areas adjacent to the endplates. 
These images were also corroborated with images from the thread tufts extended to the left and right sides. The tuft 
images showed no indication of spanwise flow on the model. Similar results were obtained for flow visualization 
performed on the flap, depicted in Fig. 8. In this case, there were small regions of spanwise flow evident near the 
endplates. This was attributed to very narrow gaps between the endplates and the flap. These narrow gaps were not 
sealed as was the main element, to allow for adjustment of the flap angle. The surface-oil flow images near the 
trailing edge of the flap were somewhat ambiguous, indicating a low-shear flow on the surface. Flow visualization 
with thread tufts did not indicate flow separation in this region. This was corroborated with boundary-layer pressure 
rake measurements, discussed in the following paragraph. Overall, the flow visualization results confirmed the 
leading-edge stall type characterized by a leading-edge laminar separation bubble and very little or no turbulent 
boundary-layer separation on the main element or flap. Furthermore, the surface flow was remarkably two- 
dimensional given the small size of the endplates and large gaps between the endplates and test-section walls. 

Flow Flow 


i i 



Figure 7. Photographs of surface-oil flow visualization the left and right sides of the model main element 
at a = 18 deg. 


Flow Flow 



Figure 8. Photographs of surface-oil flow visualization the left and right sides of the model flap at a = 18 deg. 

9 

American Institute of Aeronautics and Astronautics 


Boundary-layer velocity profiles were measured near the trailing edge of the main element and flap using the 
pressure rake shown in Fig. 3. Additional measurements were conducted with the rake positioned at y/b = ±0.25 on 
the main element trailing edge to quantify potential spanwise variation over the middle 50% of the model span. 
Data were acquired over a large angle of attack range and a subset of the velocity profiles are depicted in Fig. 9 for a 
= 8, 14 and 18 deg. These profiles illustrate the boundary-layer growth approaching stall. The profiles on the main 
element were characteristic of a turbulent boundary layer and exhibited minimal differences across the wing span. 
The velocity profiles on the flap show the classic confluence of the flap boundary-layer and the wake of the main 
element for a well-designed slotted flap system. Since the confluent boundary-layer thickness exceeded the vertical 
extent of the pressure rake, the flap velocity profiles could not be normalized by the true edge velocity. Instead the 
largest velocity measured in the profile (e.g., z* = 80 mm at a = 8 and 14 deg. and z* = 25 mm at a = 18 deg.) was 
used to normalize the profiles. The velocity profiles further support the flow visualization results that indicated an 
absence of boundary-layer separation on the main element and flap trailing edges approaching stall. 

The velocity profiles were used to calculate the integrated boundary-layer parameters described in Eq. 1-3 and 
are plotted in Fig. 10. These values could not be calculated for the flap at angles of attack greater than 8 deg. 
because the thickness of the confluent boundary layer exceeded the height of the pressure rake as described in the 
previous paragraph. The displacement and momentum thickness on the main element increased with angle of attack 
as expected, reflecting the respective increases in boundary-layer thickness and skin-friction drag. These integrated 
values exhibited more spanwise variation than was observed in the velocity profiles, but this variation is not 
considered to be significant. The boundary-layer shape factor calculated for the main element decreases with angle 
of attack thus indicating an increasing resistance to boundary-layer separation at these locations. 

The experiments with the clean, baseline model confirmed the leading-edge stall type of the airfoil including the 
presence of the laminar separation bubble and minimal or no separation on the trailing edge of the main element or 
flap. The lack of separation was further confirmed from the evolution of the boundary-layer velocity profiles with 
increasing angle of attack. Furthermore, the surface flow was found to be two-dimensional leading up to stall, 
which is an important observation given the unique installation between the small endplates. Therefore, the baseline 
aerodynamics of the model showed no anomalies that could adversely affect the evaluation of anti-icing fluids and 
ice-pellet contamination. 


© — 

— Main Element, y/b = 0.0 

A 

— Main Element, y/b = ±0.25 

B 

— Main Element, y/b = -0.25 

e — 

— Flap, y/b = 0.0 



Figure 9. Comparison of boundary-layer velocity profiles taken at the trailing edge of the main element and 
flap for three angles of attack. 


10 

American Institute of Aeronautics and Astronautics 




a (deg) 


© — 

Main Element, y/b = 0.0 

A 

— Main Element, y/b = +0.25 

B 

Main Element, y/b = -0.25 

e — 

Flap, y/b ~ 0.0 


1 A _ 

0 19- 







| 12 : 
^ 10- 







cn 

cn 

<D 

M 8 : 

O 

3 













“ o 

1 J 







<D 

a : 

O 9_ 

2 I 














0 H — i — i — i — i — i — i — i — i — i — i — i — 

-4 0 4 8 12 16 20 

a (deg) 



a (deg) 


Figure 10. Variation of integrated boundary-layer parameters with angle of attack at the trailing edge of the 
main element and flap. 


B. Effect of Simulated Leading-Edge Ice Roughness and Upper-Surface Frost 

Experiments were conducted with simulated leading-edge ice roughness and upper-surface frost in order to 
better understand the aerodynamic origins of the adverse performance effects due to the anti-icing fluids and ice- 
pellet contamination. The effect of varying roughness size on the airfoil performance is shown in Fig. 11. Bragg et 
al. 28 describe the aerodynamic effect of leading-edge ice roughness. Even though the ice roughness height tends to 
be many times larger than the local boundary-layer thickness the scale of separated flow tends to be limited to the 
immediate vicinity of the individual roughness elements. Therefore there is only a small effect on the inviscid 
flowfield until the airfoil stalls. These effects are illustrated in Fig. 11 where there was little change in the lift-curve 
slope or zero-lift angle. The chief effects of the roughness were reduced maximum lift coefficient and stalling 
angle. There was minimal change in the pitching moment characteristics prior to stall. There was no observed 
change in the drag coefficient prior to stall due to the large contribution of induced drag relative to the profile drag 
because of the unique model installation between small endplates. A significant increase in profile drag would be 
expected based upon airfoil test results. 28 The lift and pitching moment data seem to indicate that the leading-edge 
stall characteristics of the wing section were preserved even with roughness applied to the leading edge. This is 
evidenced by the sharp drop in lift after stall (particularly for the k/c = 0.00005 and 0.00010 roughness) and break in 
the pitching moment. Such behavior is not uncommon with other leading-edge stall airfoils such as the NACA 
230 12. 29 The magnitude of the reduction in maximum lift coefficient range from 20 to 30%, which is also lower 
than values reported for the NACA 23012 airfoil at similar Reynolds numbers. 29 

These results have been compiled and compared against similar tests from the icing aerodynamics literature. For 
example, Lynch and Khodadoust 30 provided an extensive review of icing aerodynamics and compiled some of the 
data shown in Fig. 12. Plotted is the effect of leading-edge roughness height on loss in maximum lift coefficient. 
Inspection of the figure shows that the maximum lift degradation measured on the thin, high-performance wing in 
the PIWT is within the range of many other airfoil studies. ’ " This comparison provides further evidence of 
suitable performance of this wing for anti-icing fluid research. 
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Figure 11. Effect of sandpaper roughness applied to wing leading edge from x/c = 0.08 on the lower surface to 
x/c = 0.05 on the upper surface. 



Figure 12. Effect of leading-edge roughness size on loss of maximum lift; comparison of present data to 
archival compilation adapted from Lynch and Khodadoust 30 with additional data from Broeren et al., 29 
Abbott and Turner 31 and Lynch et al. 32,33 


A number of studies have been conducted using roughness located on the entire wing upper surface to simulate 
ground frost conditions. 15,17,34 A comparison of this effect to roughness located on the leading edge only is shown in 
Fig. 13 for k/c = 0.00010. The main effects were observed in lift and pitching-moment coefficient. While 
significant changes in profile drag are expected, they cannot be quantified from the drag coefficient shown in Fig. 13 
because of the large contribution of induced drag and drag due to the endplates. The shift in both lift and pitching- 
moment coefficient curves for the entire upper-surface roughness coverage (i.e. x/c = 0.08 lower to x/c =1.0 upper) 
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results from the roughness-induced thickening of the boundary layer. 14,17 This increase in boundary-layer 
displacement thickness effectively de-cambered the wing which caused the reduction in lift coefficient and increase 
in pitching moment coefficient for a given angle of attack. 14,17 For this wing section, the extension of the leading- 
edge roughness (Fig. 11) over the entire upper surface (Fig. 13) did not affect the stalling angle which was about 10 
deg. for both cases. This may be due to the strong leading-edge stall tendency of this airfoil. The addition of the k/c 
= 0.00010 roughness, even over the entire upper surface, does not appear to have caused significant trailing-edge 
separation. This assumption is deduced primarily from the pitching-moment behavior, which was sharply increased 
at stall for the roughness configurations. A more gradual change in pitching moment would be expected if trailing- 
edge separation were present. 

This result was investigated further by systematically removing the upper-surface roughness beginning from the 
leading edge and working downstream. As shown in Fig. 13, the aerodynamic performance tends to be similar for 
all of the configurations once the roughness is removed from the leading edge. There appears to be a break in the 
lift and pitching-moment curves at a ~ 12 deg. for the k/c = 0.00010 roughness located on the aft 70% chord (i.e., x/c 
= 0.3 upper to x/c =1.0 upper) and aft 40% chord (i.e., x/c = 0.6 upper to x/c =1.0 upper). This break may indicate 
possible trailing-edge separation on the flap. There does not appear to be a similar break for the roughness located 
on the flap only. In all three of these cases (roughness located on aft 70% of chord, 40% of chord and flap only), the 
stalling angles are nearly identical with a small variation in maximum lift coefficient. 

The effect of the simulated frost coverage on aerodynamic performance was compared against results from 
Runyan et al. 17 who performed similar experiments on an 9.1% scale, 3-D half-plane model of a Boeing 737- 
200 AD V airplane. The experiments were conducted at near operational airspeeds, but the Reynolds number based 
upon average chord length was only about l.OxlO 6 due to the small model scale. Depicted in Fig. 14 are the lift and 
pitching-moment characteristics plotted in similar form to the original Runyan et al. 17 data. The general form of the 
performance effects of the 40%, 70% and 100% chordwise coverage is similar between the data sets. These 
similarities are more readily visible in the pitching-moment coefficient variation with lift coefficient. As described 
earlier in this section, the large chordwise extent of the roughness has the effect of increasing the boundary-layer 
displacement thickness on the aft section of the wing. The aerodynamic result is a decrease in wing camber and a 
corresponding increase in pitching moment. There is less of an effect of the roughness on lift coefficient for the 
9.1% scale B737-200ADV wing which could be due to the 3-D configuration versus the 2-D model configuration 
for the present tests. In addition, there could be Reynolds number effects for the 9.1% scale model data that would 
minimize the effect of the roughness on lift coefficient. 
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Figure 14 Comparison of aerodynamic effects of chordwise coverage of sandpaper roughness for present 
data and data from Runyan et al. 17 on 9.1% scale 3-D B737-200ADV half-plane model. 
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The simulated frost effects were compared to the effects of anti-icing fluid and ice-pellet contamination, also 
following the method of Runyan et al. 17 These results are shown in Fig. 15 for the k/c = 0.00010 roughness located 
on the flap compared to an anti-icing fluid tested at -5 °C and the same anti-icing fluid exposed to moderate ice- 
pellet precipitation at -8.4 °C. Data for the latter cases (with fluid and contamination) were obtained following the 
procedure simulating a takeoff acceleration profile and subsequent angle of attack rotation at 2.7 deg./sec, while the 
roughness configuration was taken at steady-state conditions in pitch-pause mode. The obvious effect of the 
simulated takeoff profile is that the fluid/contamination flow off of the wing during the run. This was observed 
visually and recorded during the experiments. The aerodynamic effect of this is seen in both the lift and pitch- 
moment variation with angle of attack. The linear portion of the lift curves for the two fluid cases had slightly 
higher slope than for the roughness configuration. The fluid cases caused a larger aerodynamic penalty relative to 
the clean, dry wing at lower angle of attack. As the angle of attack was increased, more fluid/contamination flowed 
off the wing thus reducing the aerodynamic impact. The pitching-moment was affected in an analogous way. For 
the fluid/contamination cases, the pitching-moment coefficient decreased more than for the roughness configuration 
as the angle of attack was increased due to the dynamics of the fluid flowing off of the wing. What is perhaps more 
interesting is the difference in stall behavior between the roughness configuration and the fluid cases. Stall for the 
latter cases occurred at lower angle of attack and lift coefficient compared to the roughness configuration. 
Furthermore, there was no clear “break” (or slope change) in the lift curve for the fluid cases. This appears to 
indicate that stall may be driven by residual fluid and/or contamination located in the leading-edge region as 
opposed to the contamination on the aft portion of the wing and flap. This is often referred to as a “secondary wave” 
of fluid and is addressed in Section III.D below. 

The data in Fig. 15 were compared once again to the Runyan et al. 17 results plotted in a common format in Fig. 
16. The lift coefficient variation with angle of attack shows similar trends between the two data sets. In the Runyan 
et al. data on the 9.1% scale B737-200ADV half-plane model, the differences in lift-curve slope between the 
transitory fluid results and the static roughness configuration are more obvious than in the present data. However, 
there are very comparable effects on maximum lift coefficient and stalling angle. The trends in pitching-moment 
coefficient versus lift coefficient are also very comparable between the two data sets. The comparisons shown in 
Figs. 14 and 16 are significant because they illustrate similar trends of fluid effects on aerodynamic performance for 
two different wing configurations. This provides further evidence suggesting that the present PIWT wing 
configuration is suitable for the fluid and contamination research. 



a (deg) 


a (deg) 


Figure 15. Comparison of aerodynamic effects of k/c = 0.00010 sandpaper roughness applied to the airfoil 
upper surface on the wing flap with fluid and fluid plus ice-pellet contamination. 
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Figure 16 Comparison of aerodynamic effects of chordwise coverage of sandpaper roughness and fluid for 
present data and data from Runyan et al. 17 on 9.1% scale 3-D B737-200ADV half-plane model. 


C. Effect of Contamination on the Flap 

A series of experiments were conducted to investigate the aerodynamic performance effects of contamination on 
the wing flap. Referring back to Fig. 15, it was observed that there was a significant difference in stalling angle and 
maximum lift coefficient between the k/c = 0.00010 roughness on the flap versus the fluid/contamination 
configurations. One explanation for this difference is that the fluid/contamination presented a more significant flow 
disturbance on the flap than the k/c = 0.00010 roughness. This sensitivity was investigated by applying much larger 
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k/c = 0.00023 roughness on the flap. The aerodynamic effect is summarized in Fig. 17. Consistent with previous 
results, the data show that the larger roughness caused a further thickening of the boundary layer and greater de- 
cambering of the wing. 14,17 This could also be interpreted as decreasing effectiveness of the flap, or as an effective 
decrease in the flap deflection. While the larger roughness also caused a significant decrease in maximum lift 
coefficient relative to the smaller roughness, the stalling angles were nearly identical at 20.4 deg. This indicates that 
the reduced stalling angle associated with the fluid/contamination cases in Fig. 15 was likely not caused by 
accumulation of fluid/contamination on the wing flap. 

The effect of the roughness on the wing flap on the boundary layer was investigated further. The boundary-layer 
pressure wake was located near the trailing edge of the main element at the midspan station (y/b = 0.0). The velocity 
profiles in Fig. 18 show that the roughness applied to the flap had a small upstream influence on the main element 
boundary layer. The effect increased with angle of attack as indicated in the profiles measured at a = 18 deg., where 
the boundary-layers for the roughness configurations were approximately 6 to 8 mm thicker than for the clean 
configuration. Larger differences in the boundary-layer profiles were expected given the large differences in 
corresponding aerodynamic performance shown in Fig. 17. 

The measurements were repeated with the pressure rake located near the trailing-edge of the flap at the midspan 
station (y/b = 0.0). These profiles, depicted in Fig. 19, illustrate the expected significant differences due to the 
roughness. The boundary layer on the flap was much thicker with the roughness meaning that the rake was not large 
enough to measure the complete velocity profile. Since the boundary-layer edge velocity was not measured, the data 
in Fig. 19 are normalized by the freestream velocity. The format of the plots is also changed to show the evolution 
of the boundary-layer profiles with angle of attack for each configuration. For the clean wing, the flap boundary- 
layer development was clearly seen for z* < 30 mm with the confluence of the main element wake for z’ > 30 mm. 
It appears that the edge of the boundary layer was measured for a = 6 deg. and that the velocity was increasing 
toward the edge for the other angles of attack for z*> 80 mm. With the k/c = 0.00010 roughness applied to the flap, 
the velocities for z’ < 30 mm are significantly reduced from the clean configuration and the profiles appear to be on 
the verge of separation near the surface. Farther out in the boundary layer the velocities appear to be increasing 
toward the edge at a = 6 and 10 deg. for z’ > 80 mm. The velocity profiles for the k/c = 0.00023 configuration are 
not significantly different than the k/c = 0.00010 profiles over the range of z* that could be measured with the 
pressure rake. Based upon the performance data it was expected that k/c = 0.00023 profiles would show a thicker 
boundary-layer than the k/c = 0.00010 profiles. Some evidence of this is shown at a = 6 and 10 deg. for z’> 50 mm. 
Taken together, these data indicate that effect of flow disturbances located on the flap, tend to manifest themselves 
primarily on the flap itself with limited upstream influence. 
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Figure 17. Aerodynamic sensitivity to roughness size on the wing flap. 
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Figure 18. Effect of roughness applied to wing flap on boundary-layer velocity profiles measured near the 
trailing edge of the main element at the midspan station (y/b = 0.0). 
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Figure 19. Effect of roughness applied to wing flap on boundary-layer velocity profiles measured near the 
trailing edge of the flap at the midspan station (p/Z> = 0.0). 
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D. Effect of Leading-Edge Flow Disturbances 

An open question that arose from the discussion of the results in Fig. 15 was the origin of the significant 
difference in stalling angle and maximum lift coefficient between the k/c = 0.00010 roughness on the flap versus the 
fluid/contamination configurations. It was concluded in Section III.C that this difference in performance was not 
due to a different level of flow disturbance on the flap. The obvious alternative explanation is that there was some 
flow disturbance near the wing leading edge for the fluid/contamination cases that caused the lower stalling angle 
and maximum lift coefficient. As described in Section III.B, the fluid/contamination was applied to the wing at a = 
-2 deg. The wing was then subjected to a simulated takeoff where the airspeed was increased to 80 knots before the 
wing was pitched through stall which is a dynamic process. 

The transitory nature of anti-icing fluid behavior suggests that some amount of fluid may remain on a wing 
surface near the leading edge through takeoff. Hill and Zierten 14 described this situation based upon their flight tests 
with a Boeing 737-200ADV airplane. The fluid applied to the wing surface naturally flows around the leading edge 
to the lower surface due to gravity and surface tension forces prior to the initial take-off acceleration. As the 
airplane accelerates, the shearing forces near the leading edge are relatively low due to the low angle of attack and 
low initial velocities during the acceleration. There is some shearing of the fluid as the airplane accelerates resulting 
in the primary wave of fluid flowing downstream. As the airplane rotates, the shear forces near the leading edge 
increase significantly and the attachment point moves farther aft on the lower surface. These effects combine to 
form the secondary wave of fluid that flows around the leading edge and downstream after rotation. While the 
secondary wave of fluid is generally much thinner than the primary wave, it can have a significant effect on 
aerodynamic performance because it is located close the wing leading edge at higher angles of attack. 

In the present series of tests, this secondary wave was simulated with 2.0 in (50.8 mm) wide layers of aluminum 
speed tape with the forward edge placed at x/c = 0.02 on the upper surface. This location was selected to be 
upstream of the laminar boundary-layer separation location that was documented with the surface oil flow 
visualization described in Section III. A. The effect of increasing tape thickness is illustrated in Fig. 20 (cf. Table 4 
for dimensional thickness range). Little or no effect on maximum lift and stall angle was observed until the height 
of the flow disturbance reached k/c = 0.00042, where there was a significant decrease in stalling angle and 
maximum lift coefficient. Incrementally increasing the disturbance height had no further adverse effect on the 
stalling characteristics up to k/c = 0.00069. The range of maximum lift coefficient was 1.91 to 1.95 at a vto// = 15.3 
deg. for the k/c = 0.00042 to 0.00069 cases. This compares to a range of C L>max = 1.91 to 1.93 and a vto// = 15.5 to 
16.0 deg. for the fluid and fluid with moderate ice pellets cases in Fig. 15. The fact that these values are nearly 
identical strongly suggests that the reductions in maximum lift coefficient and stalling angle for the 
fluid/contamination cases were at least partially due to secondary wave effects. 

Additional fluid tests were performed using a modified ramp and rotation procedure to isolate secondary wave 
effects from the overall fluid effects. Anti-icing fluid was applied to the wing using the standard procedures. The 
airspeed was accelerated to 80 knots with the wing at a = -2 deg. and then stopped to achieve fluid distribution on 
the wing that typically occurs prior to rotation. Fluid was then removed from the upper and lower surface of the 
wing aft of x/c = 0.17. In other tests, fluid was removed from the upper and lower surface separately. The airspeed 
was then accelerated to 80 knots followed by the wing rotation from a = -2 to 22 deg. Additional tests were 
performed where fluid was manually applied to the wing upper surface upstream of x/c = 0.17 (identified as “Wing 
Section A”) and to the wing lower surface upstream of x/c = 0.17 (identified as “Wing Section C”). The resulting 
aerodynamic effects are summarized in Fig. 21. The lift and pitching-moment data show only minor deviation from 
the clean wing values for angles of attack less than 12 deg. This result is expected since in all cases there was no 
fluid on the wing aft of x/c = 0.17 and therefore minimal effect on the downstream boundary layer. This result is 
similar to the static simulation tests with leading-edge roughness (Fig. 11) and aluminum speed tape layers (Fig. 20). 
Also similar is the effect on the stall angle and maximum lift coefficient that were both significantly reduced from 
the clean wing values. The lowest values of a sta u and C L max in Fig. 21 occurred when the fluid was manually applied 
to the wing leading edge on the lower surface (“Wing Section C). Fluid located in this region was likely directed 
toward the upper surface by the stagnation point moving downstream as the wing was pitched from a = -2 to 22 deg. 
This is exactly the scenario described by Hill and Zierten 14 as a secondary wave. The values of a sta u and C L>max for 
the fluid manually applied to Wing Section C configuration agree very well with the simulated secondary wave 
results presented in Fig. 20 for k/c > 0.00042. Returning to Fig. 21, the next highest values of a sta u and C L max are for 
the fluid removed aft of x/c = 0.17 after ramp. This configuration likely had less total fluid in the leading edge 
region, but fluid still remained on the lower surface after the ramp. When fluid was manually applied only to the 
upper-surface, leading-edge area (Wing Section A), the values of a sta u and C L>max were the highest. There was still a 
significant performance penalty due to fluid in this location. These fluid test results were consistent with the results 
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simulating the secondary wave. Taken together, it was concluded that the changes in stalling angle due to the fluid 
and fluid plus contamination were at least partially due to secondary wave effects. 

Another factor that may have influenced the effect of the secondary wave was the laminar separation bubble 
discussed in Section III. A. It is possible that this separation bubble could form with fluid on the surface as the wing 
is rotated through stall. The formation of the separation bubble would further promote accumulation of fluid in a 
critical location near the wing leading edge. This would increase the height of the secondary wave and contribute to 
the observed reduction in maximum lift and stalling angle for fluid/contamination cases such as in Fig. 15. This 
phenomenon appeared to occur in a small subset of the photographic documentation of the fluid tests. 
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Figure 20. Effect of flow disturbances near leading edge of wing simulating fluid secondary wave. 



Figure 21. Comparison of secondary wave effects from fluid tests. 
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IV. Summary and Conclusion 

The FAA has worked with Transport Canada and others to develop allowance times for aircraft operating in ice- 
pellet precipitation. Wind-tunnel testing has been carried out to better understand the flowoff characteristics and 
resulting aerodynamic effects of anti-icing fluids contaminated with ice pellets. The objective of this paper is to 
characterize the aerodynamic behavior of the thin, high-performance wing that is important to understanding the 
adverse aerodynamic effects of anti-icing fluids and ice-pellet contamination. All experiments were carried out at 
the NRC PIWT using the 6 ft (1.82 m) chord thin, high-performance wing section. The wing section had a span of 
7.9 ft (2.6 m) and had a single-element, slotted flap fixed at 20 deg. for these experiments. Aerodynamic 
performance data, boundary-layer surveys and flow visualization were conducted at a Reynolds number of 
approximately 6.0><10 6 and a Mach number of 0.12. 

The experiments with the clean, baseline model confirmed the leading-edge stall type of the airfoil including the 
presence of the laminar separation bubble and minimal or no separation on the trailing edge of the main element or 
flap. The lack of separation was further confirmed from the evolution of the boundary-layer velocity profiles with 
increasing angle of attack. Furthermore, the surface flow was found to be two-dimensional leading up to stall, 
which is an important observation given the unique installation between the small endplates. The evolution of the 
boundary-layer velocity profiles with angle of attack also indicated the lack of potential separation at the trailing 
edge of the wing main element and flap. 

The effect of simulated leading-edge ice roughness on maximum lift coefficient and stalling angle was large, but 
very comparable to archival data found in the technical literature. Further tests also showed that the aerodynamic 
effect of the anti-icing fluids (and fluid with ice-pellet contamination) were similar to various chordwise coverage of 
simulated frost. Experiments were also conducted with two different roughness ( k/c = 0.0001 and 0.00023) and 
sizes applied to the flap. The larger roughness size caused a slightly lower maximum lift coefficient but at the same 
stall angle as for the smaller roughness size. Boundary-layer surveys near the trailing edge of the main element 
showed that the roughness applied to the flap had a small upstream influence on the main element boundary layer. 
Overall, these aerodynamic data indicated that the primary aerodynamic effect of the anti-icing fluid (and fluid with 
contamination) was essentially a de-cambering of the wing due to increased boundary-layer displacement thickness 
on downstream portions of the wing including the flap. 

Tests were also conducted to investigate effects of a secondary fluid wave that forms as the wing is rotated in the 
simulated takeoff profile. Layers of aluminum speed tape were applied to the upper surface of the wing at x/c = 0.02 
to simulate a secondary wave of fluid. It was found that a height of k/c = 0.00042 to 0.00069 resulted in maximum 
lift coefficients ranging from 1.91 to 1.95 compared to the clean wing value of 2.2. Correspondingly, the stall angle 
was reduced to 15.3 deg. for all cases compared to the clean value of 20 deg. Tests were also conducted with anti- 
icing fluids using modified procedures to isolate the effects of the secondary wave. These tests resulted in 
maximum lift coefficients ranging from 1.94 to 2.02 and stall angles ranging from 15.5 deg. to 17.0 deg. Based 
upon the comparison of these results it was concluded that secondary wave effects could have a significant impact 
on the maximum lift coefficient and stall angle for anti-icing fluid (and fluid with contamination) tests on the thin, 
high-performance wing. 

This aerodynamic characterization of the thin, high-performance wing in the NRC PIWT has yielded important 
information about the suitability of this model for anti-icing fluid and ice-pellet contamination testing. The clean 
baseline aerodynamics of the model were consistent with expected 2-D aerodynamics and showed no anomalies that 
could adversely affect the evaluation of anti-icing fluids and ice-pellet contamination. Tests conducted with 
roughness and leading-edge flow disturbances helped to explain the aerodynamic impact of the anti-icing fluids and 
contamination. In the linear portion of the lift curve, the primary aerodynamic effect is the thickening of the 
downstream boundary layer due to the accumulation of fluid and contamination. This causes a reduction in lift 
coefficient and increase in pitching moment due to an effective de-cambering of the wing. The stalling 
characteristics of the wing with fluid and contamination appear to be driven at least partially by the effects of a 
secondary wave of fluid that forms near the leading edge as the wing is rotated in the simulated takeoff profile. 
These results have provided a much more complete understanding of the adverse aerodynamic effects of anti-icing 
fluids and ice-pellet contamination on this wing. This is important since these results are used, in part, to develop 
the ice-pellet allowance times that are applicable to many different airplanes. 
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